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An electrophoretic analysis of histone H1 of Pisum sativum L. was carried out using the collection 
of 883 accessions of cultivated peas originating from different regions of the Old World. The closely 
linked genes of five H1 subtypes (2-6) form a cluster, the gene of subtype 1 being located on 
another chromosome. The fast allelic variant of subtype 1 was not observed to the north of the 44th 
parallel. The frequency of allele 1 of subtype 5 displays a latitudinal cline and is strongly correlated 
with the sum of aerial temperatures over the vegetational period. Allele 1 of subtype 6 prevails all 
over the Old World except Central Asia and China, where allele 3 predominates. From this area 
allele 3 frequency forms a declining gradient. Alleles of subtypes 2, 3 and 4 exhibit no regularities in 
their geographic patterns. The data indicate that alleles of H1 subtype 5 and, possibly, subtype 1 in 
garden pea were subjected to climatically-dependent natural selection under conditions of primitive 
farming (without conscious selection). 
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tive trait genes. 


Introduction 


Quantitative trait genes are thought to be of primary 
importance for adaptive evolution because allelic sub- 
stitutions in them make fine alterations to the pheno- 
type (Mather & Jinks, 1982). The quantitative trait 
genes with a high degree of pleiotropy are of great 
interest as they could be involved in various adapta- 
tional processes. It is reasonable to search for such 
genes among factors influencing the higher levels of 
gene regulation. From this point of view special atten- 
tion should be paid to histone H1 genes. Histone H1 
plays a key role in chromatin folding into higher-order 
structures (Thoma et al., 1979; Allan et al., 1980), and 
is thought to be a non-specific repressor of gene 
expression (Weintraub, 1985). Besides this, histone H1 
is supposed to compete with specific gene activators 
for the regulatory sites of genes (Brown, 1984; 
Zlatanova, 1990). Recently it was proposed that 
histone H1 mediates the action of nuclear receptors as 
a participant in the induced regulatory cascade 
(Oikarinen, 1991). Therefore, alterations of molecular 
structure of histone H1 affecting its affinity to DNA or 
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its interaction with other chromosomal proteins might 
result in changes in differential gene activity. 

Histone H1, in contrast to conservative core 
histones, is highly variable in different organisms (Wells 
& McBride, 1989). Moreover, H1 is frequently repre- 
sented by several major subtypes, which may differ in 
strength of interaction with DNA (Cole, 1984) and 
differentially reduce the number of transcription initia- 
tion sites (Hannon et al., 1984). A preferential location 
of one of the H1 subtypes to certain bands of polytene 
chromosomes has been shown in Chironomus thummi 
(Mohr e¢ al., 1989). 

Our previous data suggest a possible role of H1 
genes in adaptive evolution. The rate of molecular 
evolution of histone H1 appeared to correlate with the 
rate of taxonomic evolution in insects (Berdnikov et al., 
1984). In a population of Vicia unijuga A. Br., occupy- 
ing a small territory, the allelic frequencies of one of the 
histone H1 subtypes exhibited a radial cline pre- 
sumably formed under pressure of some human- 
dependent factor (Berdnikov et al, 1992). A 
comparison of isogenic lines of garden pea that differ 
in histone H1 variants showed a mild effect of allelic 
substitutions in H1 genes over a number of quantitative 
traits (Bogdanova et al., 1989). 
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In the present work we have attempted to obtain 
further evidence of susceptibility of H1 allelic variants 
to natural selection. For this purpose we studied the 
allelic composition of H1 genes in aboriginal forms of 
garden pea Pisum sativum L. not affected by conscious 
selection. This objective was chosen for two main 
advantages. First, the pea was domesticated some 
9000-10,000 years ago (Zohary & Hopf, 1973) and in 
the course of its cultivation has occupied a wide range 
of environments over a vast territory. Second, histone 
H1 of the pea is represented by a number of electro- 
phoretic subtypes, each of them demonstrating allelic 
polymorphism that allowed genetic mapping of the 
corresponding genes. We revealed different geographic 
patterns for several H1 genes, suggesting that at least 
one of them was subjected to natural selection in the 
course of pea expansion over the Old World. 


Materials and methods 


fsolation and electrophoresis of histone H1 


Histone H1 was isolated by Johns’ (1964) method with 
further modifications (Rozov et al., 1986): About 1-2 g 
of pea leaves were homogenized in 20 ml of 0.15 M 
NaCl, the homogenate was filtered through two layers 
of gauze and centrifuged at 1500 g for 5 min. Histone 
H1 was extracted by resuspending the pellet in 5 ml of 
5 per cent HCIO,. After centrifugation the protein was 
recovered from the supernatant by adding six volumes 
of cold acetone and sulphuric acid to a concentration 
of 0.5 M. The precipitated protein was centrifuged and 
then dissolved in 0.2 ml of a medium containing 0.9 M 
acetic acid, 8 M urea, and 15 per cent sucrose. The 
preparations were subjected to electrophoresis in long 
(up to 40 cm) slabs of 15 per cent polyacrylamide/0.5 
per cent N,N’-methylenebisacrylamide gel containing 
8 Murea and 0.9 acetic acid following a modification 
(Berdnikov & Gorel, 1975) of Panyim and 
Chalkley’st(1969) method. After electrophoresis, the 
gels were stained in 0.01 per cent Coomassie R-250 in 
0.9 m acetic acid. 

The resolution of slightly differing allelic variants of 
the slowest subtype was facilitated by cleavage of the 
histone H1 molecule at tyrosine residues by treating 
with N-bromosuccinimide (1 mg ml~!) for 15 min 
(Sherod ef al., 1974). 


Histone H?/ electrophoretic spectrum 


Electrophoretic patterns of histone H1 in garden pea 
comprise, as a rule, seven subtypes (Fig. 1). The 
subtypes are denoted with numerals according to the 
increase in their electrophoretic mobility. Subtype 7 is 


detectable only in the growing parts of a plant, so it was 
excluded from our analysis. 

Identification of the proteins in question as histone 
H1 is based on extractability from chromatin with 5 
per cent perchloric acid and on the analysis of their 
molecular structure. These proteins contain an 
extremely large number of lysine residues (about 
27-30 per cent as estimated by a method of incom- 
plete succinylation) localized mostly in C-terminal part 
of the molecule that is characteristic for H1 (Belyaev & 
Berdnikov, 1985). The cleavage with N-bromo- 
succinimide and chymotrypsin (Rozov et al., 1986) 
revealed the aromatic residues (tyrosine and pheny- 
lalanine) to be located in positions characteristic for 
H1 histones (Wells & McBride, 1989). 

Each subtype can be represented by a number of 
electrophoretic variants. The variants of any subtype 
exhibit a mendelian segregation as being alleles of a 
single gene. The gene His] of subtype 1 was mapped 
on chromosome 5 (Rozov et al., 1986; Smirnova et al., 
1990); other H1 genes were on chromosone 1. The 
genes of subtypes 2-6 are linked very closely forming a 
gene cluster His(2-6) (Belyaev & Berdnikov, 1981), the 


gene of subtype 7 being located 33 cM away (A. I. 


Belyaev, personal communication). 

An allelic combination of genes of cluster His(2-6) 
attributed to a single chromosome will be referred to as 
a ‘haplotype’. An allelic composition of a haplotype is 
designated by a haplotypic formula. As about 99 per 
cent of peas have the same variant of subtype 2, its 
allelic state is not included in the formula, which is four 
symbols referring to subtypes 3-6. In accordance with 
the increase of their electrophoretic mobility the 
symbol of subtype 3 is placed first, that of subtype 4 is 
second, and so on. Each symbol is a digit correspond- 
ing to a certain allelic variant of the relevant subtype, 
the variants being numbered according to the increase 
in their electrophoretic mobility. For example, record 
1221 means that the gene of the subtype 3 (the first 
position in the formula) is represented by the allele 
encoding the slowest variant, 1, the gene of the subtype 
4 (the second position), by the allele of the second slow 
variant, 2, and so on. Sometimes some H1 subtype is 
missing; in the absence of an electrophoretic band the 
subtype behaves as null allele, so it is denoted as 0. The 
symbol X (used later in the text) implies that the 
relevant gene can be represented by any possible allele. 


Materials 


The 883 accessions of cultivated pea subspecies 
(Pisum sativum sp. sativum L., P. s. ssp. asiaticum Govo- 
rov, and P. s. ssp. transcaucasicum (Govorov) Maka- 
sheva, from the V.LR. (Vavilov All-Union Institute of 
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Fig. 1. Electrophoretic composition of pea ( Pisum sativum L.) histone H1: examples of the most frequent haplotypes. Figures 
denote H1 subtypes; superscripts, allelic variants. (Note that allelic variant 2 of subtype 3 is faster than variant 1 of subtype 4.) 15 
denotes the ‘slow’ variants of subtype 1. (a) V.LR.-2303, Arkhangelsk region; (b) V.LR.-5742, Bulgaria; (c) V.LR.-4204, Finland; 
(d) V.LR.-3716, Arkhangelsk region; (e) V.LR.-4473, Serbia; (f) V.LR.-4391, Daghestan; (g) V.LR.-4993, Daghestan; (h) V.LR.- 


4010, Greece; (i) V.LR.-6350, Afghanistan. 


Plant Breeding, Leningrad)) collection were analysed. 
For our analysis we tried to choose local forms not 
affected by systematic selection. First, we took abori- 
ginal peas collected by N. I. Vavilov and his colleagues 
(mostly before 1929) during numerous expeditions 
into regions with primitive agriculture. Also, our 
samples were supplemented with pea landraces tradi- 
tionally used by local farmers, represented in the V.LR. 
collection by accessions with precisely specified 
geographic origin. Extensively used cultivars were not 
included in the study. 

The chosen accessions were attributed to 19 regions 
of the Old World (Fig. 2}, so that each region included 
not less than 30 accessions. The majority of accessions 
was homogeneous for the H1 pattern, although some 
of them contained several patterns. Each accession was 
represented in regional samples by two plants. 

A number of accessions of Pisum sativum wild sub- 
species, namely P. s. ssp. syriacum (Boiss. et Noé) 
Berger (= P. s. ssp. humile Boiss. et Noé), and P. s. sp. 
elatius (Bieb.) Schmalh., and a wild pea species Pisum 
fulvum Sibth. et Smith, from the VIR. collection, as 
well as some other sources, were also analysed. These 
data are not included in regional samples. We also 


excluded five accessions of a very peculiar subspecies 
(often considered as a distinct species) Pisum sativum 
sp. abyssinicum (A.Br) Berger. 


Results 


Distribution of alleles of subtype 7 


The great majority of peas examined possess the same, 
‘slow’, allelic variant of the subtype 1. The ‘fast’ allelic 
variant of this subtype has been observed only in eight 
regions, belonging mainly to the subtropical climatic 
zone (Fig, 3), These are Southern Europe, Asia Minor, 
the Caucasus, Afghanistan, Tadjikistan, India, China 
and Ethiopia. The maximal occurrence of the fast allele 
is attributed to the vicinity of the 38th parallel (not 
shown); we failed to observe this variant to the north of 
the 44th parallel. It is highly improbable that this event 
is accidental (P<0.001). Besides, we found the third 
‘superslow’ allele of subtype 1 four times: in two acces- 
sions from Georgia, one accession from Spain, and ina 
very heterogenic accession V.I.R.320 of P. s. syriacum 
from Palestine. 


202 V. A. BERDNIKOV £7 AL. 


Fig. 2. An accepted division of the Old 
World into regions. n is the number of 
accessions examined. (1) Finland and 
North-Western part of the U.S.S.R. 

(n= 38); (2) North Dvina basin 

(n= 34); (3) North-Western Europe 
(n= 38); (4) Central Europe (n= 32); 
(5) Dnieper basin (7 = 31); (6) Volga 
basin (1 = 30); (7) The Urals and West- 
ern Siberia (= 34); (8} South Siberia 
and Mongolia (7 = 40); (9) Far East 
(n= 34); (10) Iberian Peninsula 

(n= 38); (11) South Europe (n= 54); 
(12) Asia Minor (= 78); (13) Cauca- 
sus (n= 43); (14) Afghanistan (7 = 48); 
(15) Tadjikistan (m = 114); (16) Western 
China (n= 51); (17) North Africa and 
Near East (1 = 33); (18) India and 
Nepal (n= 40); (19) Ethiopia (7 = 65). 


Fig. 3. Geographic distribution of the 
frequencies of the fast variant of sub- 


type 1. 


Variants of subtype 2 


Nine accessions originating from different countries 
(one from Holland, one from Russia (Perm region), one 
from Tadjikistan, one from Afghanistan, one from 
China, one from Morocco, two from Ethiopia) turned 
out to have a slightly slower electrophoretic variant 
than subtype 2, one accession from China had a variant 
with a slightly increased mobility, and all the others had 
the same common variant. 


Distribution of haplotypes of subtypes 3-6 


Tables 1 and 2 show the occurrence of different haplo- 
types of subtypes 3-6 in the regional samples 


examined. We observed a total of 33 haplotypes. Some 
of them are common in every regional sample (Table 
1): the most frequent of these are 1121, 1221, 1021 
and, to some extent, 2121, whereas others are rare 
(Table 2). For simplicity we considered the distribution 
of the allelic frequencies of each subtype indepen- 
dently, bearing in mind the constraints exerted by link- 
age disequilibrium. 


Geographic pattern of alleles of subtype 6 


Allele 1, encoding the slowest electrophoretic variant 
of subtype 6, was found to predominate in all regions 
examined except for three Asian regions, where it was 
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Table 1 Occurrence of frequent haplotypes of histone H1 fast subtypes in examined regions 
ees 


Regions 

Haplotypes 1 2 3 4 5 6 7 8 10 11 12 13 #14 #15 #16 «17 «218 «#419 

1111 13. 14 4 10 6 4 4 8 0 0 0 0 4 6 6 0 6 0 4 
1121 21 24 28 36 13 25 32 48 32 69 40 130 24 21 49 16 24 30 = 27 
1131 2 2 0 0 2 0 5 0 0 4 0 16 0 2 0 0 0 0 2 
1211 28 «618 2 10 12 14 ~= «35 8 6 6 10 0 7 2 56 2 2 2 6 
1221 26 14 48 36 24 26 33 38 26 30 70 82 41 #18 36 #418 #66 «50 161 
1231 0 0 0 0 4 6 1 4 12 12 4 4 0 0 2 0 0 4 0 
1011 0 0 14 4 2 0 0 2 0 0 3 0 0 0 2 0 0 0 6 
1021 52. 50 46 20 34 18 34 14 50 16 48 56 12 17 36 #12 #14 «28«+16 
2121 0 4 6 0 8 0 0 3 0 11 8 10 32 8 8 4 16 0 38 
1123 0 4 4 0 13 #17 «+10 = «24 6.4 8 12 4 68 145 112 0 40 0 
1133 0 0 0 0 0 4 4 4 0 0 0 0 2 31 50 15 0 0 0 
1013 0 0 0 4 0 6 6 0 0 0 1 0 0 3 2 9 0 0 0 
1023 0 0 0 0 0 0 0 0 0 0 2 0 0 9 4 11 0 0 0 
2122 0 0 0 4 0 0 0 4 0 0 0 0 32 0 28 0 0 6 0 
N 152 136 152 128 124 120 168 160 136 152 216 312 172 192 456 204 132 160 260 


Regions are numbered as in Fig. 2. 


Nis the total number of examined chromosomes, ie. the twofold number of plants, or the fourfold number of accessions. 


replaced with allele 3. Regional frequencies of allele 3 
(Fig. 4) form an evident cline with the maximum in 
Central Asia. 

As these alleles belong to certain haplotypes, we 
consider all haplotypes carrying the same allele of sub- 
type 6 as a haplotypic group. Table 3 (see also Fig. 7) 
gives proportions of the main haplotypes within XXX3 
group for regions where they are frequent; there is an 
extremely high proportion of haplotype 1123. At the 
same time, the XXX1 group has no haplotype that pre- 
dominates anywhere (Table 4). 

In the Caucasus we observed a high frequency (19.8 
per cent) of allele 2 of subtype 6, although it was 
extremely rare in other regions (Tables 1 and 2). This 
allele appears almost exclusively in haplotype 2122. 
The fourth allele of subtype 6 was seen in two acces- 
sions from Greece, one from the Minsk region, and one 
from the Irkutsk region. 


Geographic pattern of alleles of subtype 5 


We distinguished three electrophoretic variants of sub- 
type 5. (Variant 2 is represented by one of two forms 
with hardly detectable differences in electrophoretic 
mobility: less than 0.5 per cent.) Figure 5 demonstrates 
the geographic distribution of their frequencies. Note 
that there is a clear trend for allele 1 frequency to 
increase in the northern direction. The trend is dis- 
rupted only in Tadjikistan because of the high fre- 


quency of haplotype 1211. We found that more than 
half of the specimens having this haplotype were attri- 
buted to a narrow region of Gorny Badakhshan. This 
alpine region of the Pamirs is characterized by a severe 
climate, similar to that of the northern boundary of 
agriculture. 

The similarity is well reflected by the sum of aerial 
temperatures over the vegetational period with 
temperatures above 10°C (Agricultural Atlas of the 
World, 1972). In Fig. 6a frequencies of allele 1 of sub- 
type 5 are plotted against this temperature sum (the 
sample from Tadjikistan being replaced by the sample 
from Gorny Badakhshan). This frequency decreases 
rapidly with the increase in the sum of vegetational 
temperatures. A logarithmic transformation of the axis 
of temperature sums linearizes this relationship (Fig. 
6b), the correlation coefficient being -—0.84 
(P<0.005). It should be noted that the trend becomes 
more explicit (correlation coefficient — 0.91) if we con- 
sider the proportion of XX11 haplotypes within 
haplotypic group XXX1. 


Alleles of subtypes 3 and 4 


There are two points of interest concerning these sub- 
types. First, no pea plant was found to possess the 
allelic combination 22XX. Second, the null allele of 
subtype 4 is often observed, while the null alleles of 
other subtypes are scarce. Figure 7 shows the 
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geographical distribution of the frequencies of classes 
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Fig. 4. Geographic distribution of the 
frequencies of allele 3 of H1 subtype 6. 


Table 3 Proportion of histone H1 haplotypes within XXX3 haplotypic group for samples from Central Asian regions and other 
regions combined 


Haplotypes (%) 
Samples 1123) 1133 1013 1023 1223 1213 2123 1113 1233 2133 0103 nt Nt 
Afghanistan 58.1 265 2.6 77 0.0 0.0 0.0 1.7 0.0 3.4 0.0 117. =—192 
Tadjikistan 63.3 218 09 1.7 3.5 2.6 2.6 2.6 0.8 0.0 0.0 229 456 
China 73.3 99 5.9 7.2 0.0 2.6 0.0 0.7 0.0 0.0 0.0 152 204 
Other regions combined 65.5 6.3 76 0.9 8.1 5.4 3.6 0.0 1.8 0.0 0.9 223 2678 
Total sample 65.3. 15.2 4.3 3.6 3.6 3.1 1.9 1.2 0.8 0.5 0.3 721 3530 


tnis the number of chromosomes with XXX3 haplotypes. 
+N is the total number of chromosomes (the fourfold number of accessions) examined. 


Table 4 Proportion of the most frequent histone H1 haplotypes within XXX 1 haplotypic group for combined samples 
corresponding to major parts of the Old World 


Haplotypes (%) 

Samples 1221 1021 1121 2121 1211 1011 1111 1231 1131 vt Nt 

West Europe! 30.9 21.8 29.1 4.2 47 3.5 2.4 2.7 0.7 595 648 
East Europe* 199 304 1846 19° 173 0.3 6.6 1.8 1.8 618 700 
North-East Asia? 25.55 25.5 31.9 1.2 56 08 3.2 6.4 0.0 251 296 
Central Asia‘* 224 20.2 26.7 62 186 0.6 3.7 0.6 0.6 322 852 
India® 43.9 246 26.3 0.0 18 0.0 0.0 3.5 0.0 114 160 
Anterior Asia and North Africae 33.8 164 318 108 16 0.0 1.8 0.7 2.9 560 =616 
Ethiopia’ 61.9 62 104 146 2.3 2.3 1.5 0.0 0.8 260 260 
Total sample 31.1 21.1 254 5.8 83 12 3.3 2.0 1.3 2731 3530 


Regions are numbered as in Fig. 2. 

1 Region no. 3, 4, 10, 11;? region no. 1, 2, 5, 6, 7; 3 region no. 8, 9; * region no. 14, 15, 16; ° region no. 18; ° region no. 12, 13, 17; 
7 region no. 19. 

+n is the number of chromosomes with XXX1 haplotypes. 

+Nis the total number of chromosomes (the fourfold number of accessions) examined. 
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Fig. 5. Geographic distribution of the 
frequencies of allele 1 of H1 subtype 5. 


Second, evolution of cultivated plant species under 
the pressure of natural selection should be accelerated 
compared with that of wild counterparts because of the 
more homogenous conditions provided by syn- 
chronous sowing and other techniques, leading to a 
decrease of environmental components in the pheno- 
typic variance and hence to the better exposure for the 
selective advantages of certain genotypes. Moreover, 
population sizes of cultivated plants are large enough 
for the effect of genetic drift to be neglected. For 
example, the seed stock needed for sowing the field 
typical for Highland Pamirs (Vavilov & Bukinitsh, 
1959) should be not less than 10°. 

Third, there is often archaeological and historical 
information available concerning the region, and the 
time of domestication and history of a cultivated plant. 
Finally, cultivated species have often been genetically 
well studied. 


Evolution of cultivated pea and its H1 composition 


A complex of archaeological, botanical and genetic 
data shows that the pea, together with barley, wheats 
and lentil, was among the first domesticated plants. 
The domestication event, coinciding with the ‘Neolithic 
transition’ to productive farming, occurred approxi- 
mately 10,000 years ago presumably in the region of 
the so-called Tauro-Zagros Arc (Ammerman & 
Cavalli-Sforza, 1984). Zohary and Hopf (Zohary & 
Hopf, 1973; Hopf, 1986) assume the Anatolian sub- 
species P.s. syriacum to be an ancestor of the cultivated 
pea, P.s. sativum. Radial migration of Neolithic farmers 


and their culture from this centre with a mean rate of 
approximately 1 km per year (Ammerman & Cavalli- 
Sforza, 1984), was followed by cultivated plants 
expanding over a constantly widening area. This migra- 
tion was followed by genetic differentiation of culti- 
vated species. The process lasted for several thousand 
years of farming, during which time the agrotechnical 
level was extremely low and the exchange of seeds 
between distant regions was probably non-significant. 

Some regional forms of pea, formed in the course of 
its expansion over the Old World, achieved a sub- 
species taxonomic rank (Govorov, 1937; Makasheva, 
1979). Subspecies P.s. asiaticum, characterized by 
early ripening and resistance to late frost, was isolated 
from Asia, except its Boreal part, and probably was 
formed in the Central Asian highlands. Subspecies P.s. 
transcaucasicum, a peculiar form with narrow leaflets 
and very small seeds, is common in the Caucasus, 
where it is used as a fodder crop. 

The ancestral genetic pool of the domesticated pea 
probably contained haplotypes 1121, i221, 1021 and 
2121, which still predominate in peas all over the 
world except Central Asia and China (Table 4). Two of 
these haplotypes, 1121 and 1021, were found in the 
examined specimens of P.s. syriacum - the putative 
ancestor of the cultivated pea. 

It should be noted here that an allele of any H1 gene, 
once arisen by mutation, can hardly be transferred to 
other haplotypes by cross-over, because of the close 
linkage between genes of H1 subtypes 2-6, and 
because the pea is a self-pollinator. (Rather common 
cross-pollination is reported only for peas from Afghani- 
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Fig. 6. Frequencies (p} of allele 1 of H1 subtype 5 plotted 
against the sum of aerial temperatures above + 10° (S) 
(Agroclimatic Atlas of the World, 1972) for examined 
regions. (a) Non-transformed scales, (b) scale of S is logarith- 
mically transformed. Numbers correspond to those of Fig. 2. 
The sample of Tadjikistan was replaced with that of Gorny 
Badakhshan (no. 15a) when calculating the correlation 
coefficients, 


stan (Govorov, 1928). We also found that the substan- 
tial part of Afghani and Tadjik specimens still display 
heterozygosity for histone H1 haplotypes.) 

In this respect a possible role of repeated mutations 
should not be neglected. Histone H1 genes contain a 
series of imperfect repeats encoding short amino acid 
sequences of five to six residues (Von Holt et al., 1979; 
Wells & McBride, 1989). They could promote intra- 
genic deletions and duplications as a result of slippage 
mispairing (Dover, 1989), or unequal crossover (Smith, 
1976). Thus, histone H1 genes might have an increased 
rate of mutations that alters the length of the H1 mole- 
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cule for a fixed segment. Such length alterations result- 
ing in changes of electrophoretic mobility are possibly 
responsible for repeated appearances of the same 
allelic variants. A scheme of mutational transitions can 
be suggested for the four above mentioned haplotypes: 


1021 


——~ 


2121 1121 «——_——_» 1221 


The null allele of subtype 4 can result from a loss-of- 
function mutation. 

Among three Central Asian regions rich in XXX3 
haplotypes, we found an extremely large proportion of 
haplotype 1123 (Table 3, Fig. 7). This can be easily 
explained by founder effect if one assumes the muta- 
tional transition of haplotype 1121 into 1123 to give 
rise to all the diversity of XXX3 group, where haplo- 
typic frequencies are still far from equilibrium. This 
mutation was probably fixed in some isolated popula- 
tion, which might be ancestral for subspecies asiaticum 
distributed in this region. Most likely it happened when 
a primitive farmer carried the pea to the territory of 
Afghanistan. The derivatives of this population subse- 
quently occupied an extended region of mountainous 
Central Asia and later penetrated into China which 
was devoid of this culture until the first century B.c. 
(Makasheva, 1979). Moderate latitudes from the 
Dnieper to the Yenissei were settled both by the 
peoples of Europe and of Central Asia. These migra- 
tional flows in their interaction could form the 
observed smooth gradient of XXX3 haplotypes declin- 
ing in the western direction. (The tendency of the 
Asiatic subspecies to cross-pollination can account for 
the penetration of these haplotypes in other pea sub- 
species.) 

It is intriguing that the distribution of another culti- 
vated pea subspecies, P.s. transcaucasicum, coincides 
with a high concentration of haplotype 2122. Notably, 
we found this haplotype in one of three examined 
Transcaucasian accessions of wild subspecies P.s. 
elatius. This wild pea could hardly have been an 
ancestor of the cultivated P.s. transcaucasicum. It 
seems more likely that the latter originated from some 
form of P.s. sativum, and haplotype 2122 arose from 
the 2121 by mutation. Later it could penetrate into a 
local P.s. elatius population. 

Thus, we face a coincidence of the territorial distri- 
bution of cultivated pea subspecies with the high fre- 
quencies of alleles of H1 subtype 6. Of course, there is 
no evidence for histone H1 genes participating in the 
taxonomic differentiation of pea. Nevertheless, it seems 
not to be impossible that changes of subtype 6 could 
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influence some traits undergoing selection during the 
formation of new subspecies. 


Latitudinal variation in the allelic patterns of histone 
H1 genes 


The observed correlation of the frequency of allele 1 of 
subtype 5 with the sum over the temperatures over the 
vegetational period is the most remarkable result of 
our study, as it can hardly be explained without the 
assumption of some climatically-dependent selective 
factor which is able to distinguish between the allelic 
variants of histone H1. It should be noted that only one 
of five closely linked genes exhibits a correlation 
between the allelic frequency and the climatic factor. It 
is difficult to suppose that any non-selective factor such 
as founder effect with further migration could be 
responsible for the observed pattern. Any event that 
affected the territorial distribution of pea genotypes 
most probably should change allelic frequencies of 
closely linked genes in a concerted manner. 

The observed picture seems also to exclude a hitch- 
hiking explanation, implying that the selection acted in 
fact on some unknown locus linked to the H1 haplo- 
type. In such a case the high frequency of allele 1 of 
subtype 5 in the regions with a rigorous climate would be 
from some certain haplotype of the XX1X class, which 
happened to be linked with a hypothetical advan- 
tageous gene. However, this is not the case. In two nor- 
thern regions the considered allele is found in two 
haplotypes, 1111 and 1211, with comparable frequen- 
cies, while in the Pamirs it appears in five haplotypes: 
1111, 1011, 1013, 1113 and 1211, the latter one sub- 
stantially predominating. 


\, 
Dn 
p 


Fig. 7. Geographic distribution of the 
following haplotypic classes: (1) 21X1; 
(2} 11X1; (3) 10X1; (4) 12X1; (5) 
21X3; (6) 11X3; (7) 10X3; (8) 12X3. 
The length of columns corresponds to 
the frequencies of classes in regional 
samples. 
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The northern boundary of distribution of the fast 
allelic variant of subtype 1 appeared to coincide 
roughly with that of the natural habitat of wild forms of 
pea. This might indicate that to the north of the wild 
pea area some natural selection acted against this allele, 
while polymorphism of subtype 1 in subtropical 
environmental conditions seems merely to be tolerable. 
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